A Relevant Two Qubit Bell Inequality Inequivalent to the CHSH Inequality 
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We computationally investigate the complete polytope of Bell inequalities for 2 particles with 
small numbers of possible measurements and outcomes. Our approach is limited by Pitowsky's 
connection of this problem to the computationally hard NP problem. Despite this, we find that 
there are very few relevant inequivalent inequalities for small numbers. For example, in the case 
with 3 possible 2-outcome measurements on each particle, there is just one new inequality. We 
describe mixed 2-qubit states which violate this inequality but not the CHSH. The new inequality 
also illustrates a sharing of bi-partite non-locality between three qubits: something not seen using 
the CHSH inequality. It also inspires us to discover a class of Bell inequalities with m possible 
n-outcome measurements on each particle. 



How are actions and events in different places con- 
nected to one another? Normally we imagine that the 
correlations were arranged in the past. Both my socks 
are black since I put on a pair this morning. However our 
quantum mechanical theory of the world is more compli- 
cated 1] . Correlations are created in at least one more 
way. Some possibilities are that a) correlations are ar- 
ranged through faster than light influences in the present, 
b) correlations are arranged in a many world scenario, c) 
correlations just occur: they constitute a primary con- 
cept, preventing us from consistently thinking about lo- 
cal subsystems. Since these possibilities are disliked by 
many physicists, we have studied the set of correlations 
which can be generated by the past, and when quan- 
tum mechanics goes beyond this. Our goal is first to 
find a simple set of conditions - generalised Bell inequal- 
ities - which describe the boundaries of the set of past- 
generated correlations. This set is often called the set of 
common cause correlations, local hidden variable (lhv) 
correlations, local variable correlations, or local realistic 
correlations. Our second goal is to see when quantum 
mechanics goes outside this set. 

One reason for studying this boundary so closely is the 
fundamental question: which correlations can be gener- 
ated in this way? A second reason is that a violation 
of Bell inequalities gives a signature for useful entangle- 
ment. For instance, violation of a certain Bell inequality 
by an N qubit state implies that the state is distillable U : 
perfect bipartite entanglement can be extracted from it. 
Also, Bell inequalities can be used as a simple test for the 
security of quantum cryptography 3] . The rough idea is 
that past-generated correlations could have been created 
(and thus known) by an eavesdropper, and so are not 
useful for cryptography, whereas other kinds of correla- 
tions cannot be created by the eavesdropper and so are 
useful. The connection holds very closely in some of the 
main cases of interest . 

Bell inequalities are also related to classical communi- 
cation complexity p|: how much communication do two 
parties need in order to perform some joint task? A fi- 
nal reason is that no experiment has definitively demon- 
strated correlations outside the past-generated set[f| . We 



hope to one day find an inequality which will allow us to 
do this without waiting for improved technology. 

A typical experiment to test for correlations begins by 
creating a particular quantum mechanical state of two 
particles, and sending one to site A, and the other to 
site B. We then perform one of a certain number, tua 
and tub say, of possible measurements, %a and is, at 
each of the two sites, A and B. Each measurement has 
a certain number, ua and n^, of possible outcomes, ja 
and j'b. We then repeat the experiment many times to 
get accurate probabilities for each set of joint outcomes, 
P{jA,3B\iA> The d = mAmBriAnB probabilities can 
be thought of as a point in a d dimensional space. We are 
interested in the set of all points which can be described 
using past-generated correlations. This set is convex, and 
the boundary is defined by hyperplanes. ft is straight- 
forward to list the n2 A n^ B vertices. We would like to 
know the faces, otherwise known as the Bell inequalities. 
We want to know how many different types of faces there 
are, and moreover, which ones are relevant for quantum 
mechanics. 

Characterising the set of past-generated correlations 
is difficult. Here difficult is meant in a technical sense. 
Suppose we are given a point and asked if it is in the 
set. Finding the answer is a finite computation, but 
will take more time with more possible measurements, 
niAiTiB- In fact, Pitowsky has shown this problem to be 
NP-complete [7j. Furthermore, suppose we are given an 
inequality, and wish to know whether or not it is a face. 
This problem is of similar difficulty: co-NP complete . 

Since the general problem is so hard, we calculated 
the Bell inequalities for various small numbers of mea- 
surements and outcomes. Surprisingly, we have found 
that for small numbers of measurements and outcomes 
there are very few inequivalent Bell inequalities. For the 
case (rriA = 2. tub =2,ua = 2, ns = 2) (hereafter called 
2222) , Fine [9j has shown that up to certain equivalences 
which we shall describe below, the only Bell inequality 
is the CHSH [ljj. F° r the case 3322 we have found that 
there is only a single new inequality. Considering the 
complexity of the problem, this is very surprising and 
simple. Furthermore, the inequality is relevant, since 
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there are states which violate it but do not violate the 
CHSH inequality. We believe this is the first time that 
performing more than two measurements has been shown 
to be useful for detecting non-local correlations. 

The confidence gained from such results (we shall de- 
scribe more) has also helped us to generalise this in- 
equality to the case mmnn. Thus the computational 
approach, whilst limited, has proved rather useful. How- 
ever we do not know if these inequalities are relevant: are 
there states which do not violate the previous inequali- 
ties which violate these new ones? Computationally, it 
seems the answer is "no", though we have no system- 
atic method for checking this. But if these inequalities 
are not relevant, this would further reduce the number 
of incquivalent, relevant Bell inequalities. 

Before describing our results, we must describe more 
precisely the set of correlations, and what we mean 
by equivalent and relevant inequalities. As previously 
stated, the correlations live in a d dimensional space 
with components representing P{jAi Jb|£.A) *b)- Some of 
these components are redundant, however. For exam- 
ple, for any fixed measurement, there has to be at least 
one outcome, ie. Ylj A j B P (J-AjJb|*Aj*b) = 1- Also the 
correlations we are interested in, those of quantum (and 
classical) mechanics, do not allow faster than light sig- 
nalling. In other words, the distribution of outcomes on 
one particle does not depend upon the choice of mea- 
surement made on the other: Ylj A P {3A,3B\iAi is) = 
Y]j A P{jA, Jb|^4, is)- Therefore we work in the subspace 
of dimension 

d 2 = m A m B (jiA-l){n B - l) + m A (nA-l)+m B (n B - I) 

(1) 

which satisfies all of these constraints. Our subspace 
can be labelled using the components P(jA,jB\iA, *b) 
for j A = 0..(n A - 2), j B = 0..(n B - 2), i A = l..m A , i B = 
l..m B ; the components P{jA\iA) for ja — 0..(ua — 2), 
iA = l-iriA', and the components P(j B \i B ) for j B = 
0..(n B - 2), i B = l..m B . 

We are interested in the faces of the convex set in this 
reduced space. These are "tight" inequalities. There are 
of course other inequalities which are satisfied by all the 
points inside the past-generated set, but such inequalities 
are less useful for detecting non-local correlations. An 
inequality will describe a d 2 — 1 dimensional hyperplane, 
and so has d 2 components. 

Some of these faces are equivalent. For example, when 
defining the experiment, we have to decide which out- 
come j'a or j B is which, and which measurement iA or 
i B is which, and which particle is A and which B. Since 
these choices are arbitrary, we shall consider two inequal- 
ities to be equivalent if they can be converted into one 
another simply by relabelling these local choices. 

Having found the incquivalent faces, we would like to 
know if they are violated by quantum mechanics. As 
our motivation for studying bell inequalities comes from 
quantum physics, we are interested only in the faces 
which are violated. 



Given two inequalities which are violated, we define 
the first to be non-redundant if quantum mechanics gives 
a point P(jA, jB\iA,iB) which violates the first inequal- 
ity , but which does not violate the second inequality (or 
any inequality equivalent to the second). Given several 
inequalities, we could look for the minimal set of inequal- 
ities such that none are redundant. 

We are more interested in a classification which comes 
from quantum mechanical states. Given two inequali- 
ties we define the first to be relevant if there exists a 
quantum state which violates it for some choice of mea- 
surements, but does not violate the second inequality for 
any choice of measurements. Similarly for the second in- 
equality. Given a set of inequalities, we want to find the 
minimal set of relevant ones. Note that this will lead to 
a minimal set no larger that the set of non-redundant 
inequalities: in fact it may be much smaller. 

One is often interested in quantum systems of a certain 
size, like 2-qubit systems. How many relevant inequali- 
ties are there for 2-qubits? Before the present work, only 
one - the CHSH - was known to be relevant. An open 
question was whether more measurements (or more out- 
comes) could help. Here we show that inequalities with 

3 measurements are relevant. We do not know whether 

4 or more will help. 

More outcomes may also be useful even on qubits, since 
we could perform a several outcome POVM. Whilst the 
usefulness remains an open question, we can at least put 
an upper bound of d 2 outcomes for any measurement on a 
d-dimensional system. Thus for 2-qubits, it is not useful 
to have more than 4 outcomes for any one measurement. 
The reason for this is that a POVM in a <i 2 -dimensional 
space (the dimension of the density matrix) can always 
be viewed as a classical probabilistic mixture of POVM's 
with d 2 outcomes [Tl| . In other words, the many out- 
come POVM can be viewed in two stages. The first rolls 
an independent dice to decide which few-outcome mea- 
surement to make. The second performs the few out- 
come measurement, and gives the appropriate outcome. 
Adding local randomness which is under the control of 
the lhv model cannot add non-locality, and so any non- 
locality present in such a many outcome measurement 
must already be there in the few outcome measurements, 
and hence in a few outcome inequality. 

The limit on the number of useful outcomes is also in- 
teresting since it suggests there are inequalities which are 
irrelevant for 2-qubits, but which are useful for higher 
dimensional systems. For example the 5-dimensional 
CGLMP inequality [H, which deals with the case 2255, 
is irrelevant for qubits (it has too many outcomes). On 
the other hand, there are 5-dimensional states which vi- 
olate this inequality which are not known to violate any 
lower dimensional inequalities. With our present knowl- 
edge, this inequality is indeed useful for such systems. 

In order to find all the inequivalent inequalities, we 
have several software tools. The main one is a linear 
programming tool which takes a list of vertices as in- 
put and, after some time, outputs all the faces [l3j . 
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We have written a small matlab program which, given 
rriA, Trig, riA,nB, produces a list of the vertices. The ver- 
tices are given by distributions which factor into two local 
probability distributions, ie. 



P(jA,jB\iA,iB) = P{jA\iA)P{jB\iB) 



(2) 



and for which all the local probabilities are either or 1, 
eg. for 2222, 



P(JA 


= 0\i A = 


1) 


= o, 


(3) 


P(JA 


= 0\i A = 


2) 


= 1, 


(4) 


P(JB 


= 0\t B = 


1) 


= 1, 


(5) 


P(3B 


= 0\t B = 


2) 


= 0. 


(6) 



After we have the list of faces, we put this into a second 
matlab program which removes the equivalences, leaving 
us with the inequivalent inequalities. These software are 
all deterministic, and so give the exact solution. The 
bottleneck is the freely available linear programming tool, 
which is optimized for certain kinds of convex sets, but 
not for the equivalences which we have here. We have 
a final piece of software which, given an inequality 
and either the size of the quantum system (eg. 2 qubits) 
or a specific quantum state, probabilistically finds the 
maximum value of the inequality. 

For the case 2222 the software reproduces Fine's re- 
sult that there are only two types of inequality. One 
is the trivial one that probabilities are positive, ie. 
P{.3Ai3B\iAi > 0. This occurs TO^m^n^s = 16 
times (to cover all the joint probabilities). There is no 
need for inequalities stating that probabilities should be 
no greater than 1, since this follows from all the proba- 
bilities being positive. The other type of inequality is the 
CHSH, which occurs 8 times. We write it here in a form 
closer to that of the CH inequality |l5jl 



Ichsh = PiAxBx) 
-(P(Ai) 



P{A 2 B 1 ) + P(A 1 B 2 ) - P{A 2 B 2 ) 
P{Bi)), (7) 



where P(AB) is the probability that when A and B are 
measured we get the outcome for both measurements. 
Ichsh < for lhv correlations. Quantum mechanics can 
attain values up to — | . 

It will be useful for later on to write this inequality in 
the following way: 







-i 





Ichsh — 




i 








i 





(8) 



where the table gives the coefficients we are to put in 
front of the probabilities: 



/ 



P(Bi) 
\P(B 2 ) 



PiAtBt) P{A 2 B{) 
P(A 1 B 2 ) P(A 2 B 2 ), 



(9) 



Next we computed the case 2322. This is a choice 
between two 2-outcome measurements on one particle, 
and between three 2-outcome measurements on the other 
particle. Here we found no new inequalities. We have 
only that the probabilities must be positive, and CHSH 
inequalities where results of one the 3 measurements is 
ignored, eg. 



T' 

1 CHSH 



P(A!B 2 ) 
-{P{M)- 



P(A 2 B 2 ) 
P{B 2 )). 



P{A 1 B 5 ) 



P(A 2 B 3 ) 
(10) 



There are m^msn^nB 



and 



24 " positive probability" faces, 
= 24 CHSH faces - we have to choose 2 of 



the three possible measurements for B, and once these 
are chosen we have the 8 versions of the CHSH inequality 
which appear in the 2222 case. This gives a total of 48 
faces. 

We have analytically extended this result to the case 
2to22. We find that there are no new inequalities for 
this case. Our proof is essentially to note that the proof 
of Fine for the 2222 case extends naturally to the 
2m22 case. Fine's proof works by starting with the mea- 
sured probabilities P(jA, JbKa,zb), which are assumed 
to satisfy the CHSH inequalities. He then constructs a 
lhv model which reproduces the measured probabilities. 
We shall describe the construction for the 3222 case: the 
general to222 case follows very naturally. 

First we define (3 to be the minimum of 8 quantities: 
P(-Bi), P(AiSi) + P(B 2 ) - P{A l B 2 ), and the other 6 
quantities which come from exchanging A\ for A 2 or 
A% and B\ for B 2 in the previous expressions. We 
set P(Bi,B 2 ) = (3. P(Bi) and P(B 2 ) are experimen- 
tally measurable so we can complete the distribution for 
B x and B 2 by P(B U B 2 ) = P(B{) - 0, P(fli,B 2 ) = 
P(B 2 ) - f3_, and P{B U B 2 ) = 1 - P{B 1 ) - P{B 2 ) + [3, 
where P(B) = P(B = 1). One can check that all these 
probabilities are positive, using the fact that all the mea- 
sured probabilities P{]a, Jb\ia, *b) are positive. 

We extend this to a lhv model for A\, B\ and B 2 . We 
define a to be the minimum of P(Ai, Bi), P{A\, B 2 ), 
[3 and S3 - {P{A X ) + P{B X ) + P{B 2 ) - P{A X ,B X ) - 
P(A 1: B 2 ) - 1). We set P{Ai,Bi,B 2 ) = a. We can 
check this is well defined using the CHSH inequalities. 
We complete the distribution for (A\,B\,B 2 ) using the 
quantities we already have. ie. 

P(A U B^B^ = P^Bj-a, (11) 
P{A U B 1 ,B 2 ) = P(A u B 2 )-a, (12) 
P{A U B U B 2 ) = P(A 1 )-P(Ai,fl 1 ) 

-P(A 1 ,B 2 )+a, (13) 
P(A 1 ,B 1 ,B 2 ) = p-a, (14) 
P(A 1 ,B 1 ,B 2 ) = P(B 1 )-P(A 1 ,B 1 )-{p-a), (15) 
P(A 1 ,B 1 ,B 2 ) = P{B 2 )-P{A x ,B 2 )-{p-a), (16) 
P{A\,Bi.B 2 ) = P(A 1 ,B 1 )+P(A 1 ,B 2 ) + (p-a) 

+l-P{A l )-P(B l )-P{B 2 ).(17) 
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That these are all positive follows from the CHSH in- 
equalities. 

In a similar way, we make lhv models for the triple 
{A 2 ,B X ,B 2 ), and the triple {A 3l B 1 ,B 2 ). 

We finally extend this to a distribution for 
{A 1 ,A 2 ,A 3 ,B 1 ,B 2 ) by 



P(A 1 ,A 2 ,A 3 ,B 1 ,B 2 ) 



P{A 1 \B 1 ,B 2 )P{A 2 \B 1 ,B 2 ) 
*P(A 3 \B 1 ,B 2 )P(B 1 ,B 2 ){18) 



where P(Ai|Bi,B a ) = P{A 1 ,B U B 2 )/P{B 1 ,B 2 ). This 
gives us a well defined lhv distribution which reproduces 
all the measured probabilities. 

For three possible measurements on each side, the 
case 3322, Garg and Mermin 0] have shown that the 
CHSH inequalities are not the only faces of the clas- 
sical polytope. They found a point which satisfies all 
the CHSH inequalities, but does not admit a lhv model. 
A complete list of the faces have been computed by 
Pitowsky and Svozil 0. They found 684 faces. Re- 
moving equivalent inequalities, we find a single new in- 
equality. The 684 faces of the 3322 polytope are made 
up as Tfijijiij^Tijiji r — 36 "positive probability" faces, 

8 I I =72 CHSH faces, and 576 equivalent new faces. 
The new face is 



•£3322 = 



Vo 



1 \ 



1 1 I 
1 1 -1 

1-10/ 



(19) 



This expression satisfies 13322 < for past-generated 
correlations. For quantum mechanics a numerical opti- 
mization suggests that the maximum value is 0.25. This 
value can be attained by the maximally entangled state 



^=(|o,i>-|i,o». 



(20) 



The measurements all lie in a plane, so we denote their 
position by a single angle, the angle which they make 



with the z-axis in the Bloch sphere. 
A 3 =2f,B x = 4?, B 9 = vr, and B 3 



A\ = 0, A 2 = 5, 



13 = -, -di = 3p, B 2 = 7T, and B 3 = ?f. 

The most interesting feature of this inequality is that 
there exist states which violate it which do not violate 
the CHSH inequality. For example, consider the 2-qubit 
state 



~ 0.0129. The measurements for this violation are Von- 
Neumann measurements in the directions {0 az i m ,9 po i ar ), 
where 6 a zim is the azimuthal angle with the z-axis, and 
Opoiar is the polar angle in the x-y plane, when we set the 
z-axis to be in the direction of |0)), and the x-axis to be 
in the direction of |0) + |1): 

(23) 
(24) 
(25) 
(26) 
(27) 
(28) 



Ax 


= fa,0), 


A 2 


= (w-t/,0), 


A 3 


= (0,0), 


Bi 


= (ir-X.0), 


B 2 


= (x,o), 


B 3 


= (tt.O), 



where cos 77 = -Aj= , and cos x - 

To compare the inequalities CHSH and 13322, we nu- 
merically calculated the maximum violation, Tr(Bp), of 
-£3322 for all possible Von-Neumann measurements for a 
family of states parameterised by 9, 

P9 = ^CHSHPcos8\0,0)+sin8\l.l) + (1 — ^CHSh)P\0,1) j 

(29) 

where Xchsh is chosen so that each state in the family 
pe gives the maximal value of the CHSH inequality which 
can be obtained by lhv theories. In order to give some 
meaning to the size of the violation, we have re-scaled 
Ichsh and -Z3322 so that the lhv maximum is 1, and 
the maximally mixed state p = | gives the value 0, ie. 

Ichsh = 2I C hsh + 1, ^3322 = ^3322 + 1. The results are 
in FIG 1. 




a = O.85P| 0) +0.15P| 0)1> , 
where P^ is the projector onto the state 



(21) 



and 



FIG. 1: Maximum value of two Bell inequalities for a family 
of states. The straight horizontal line is that of Ichsh, whilst 
the curve is for 73322 . 



^(2,0,0) 



|M». 



(22) 



One can check (using the Horodecki criterion [T||) that 
this state does not violate the CHSH inequality. How- 
ever it does violate the 3322 inequality, giving a value 



The software which calculates these maximum quan- 
tum mechanical values converges quickly to very consis- 
tent results, giving us confidence that they are correct. 
We see that the new inequality is most important not for 
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states near the maximally entangled state (6* = f ), but 
rather for states with less symmetry. 

It also seems that the CHSH inequality is still rele- 
vant: there are states which violate it which do not vi- 
olate i3322- This is an illusion, due to the fact that we 
only maximized the violation over non-degenerate Von- 
Neumann measurements. Surprisingly, the maximum vi- 
olation of the new inequality is often given by degenerate 
measurements. For example, if we take 13322 and set 
A3 = 1 and B\ = 1, the remaining measurements give 
us the CHSH inequality. Thus given 13322, the CHSH 
inequality is no longer relevant. 

^3322 shows a very direct non-locality in the states of 
equation l|29[) . It is worth noting that such states are 
also non-local by Popescu's "hidden non- locality" crite- 
rion [l^, [23|. In this one first makes local filtrations to 
the state on both particles, and then performs a stan- 
dard CHSH test on the state which will emerge if both 
particles pass the filters. One only looks at the data in 
the case where both filters are passed, and if this data 
violates the CHSH inequality we are assured that the 
original state is non-local. For our states we would apply 
a local filter to particle A which lets state |1) pass, and 
absorbs state |0) with high probability. We would simul- 
taneously apply a local filter to particle B which lets |1) 
pass, and absorbs |0) with high probability. The idea is 
that each component of the entangled state is only fil- 
tered once, whereas the noise term is filtered twice. If 
both particles pass the filter, the state is very close to a 
pure entangled state. Since all pure qubit states violate 
the CHSH inequality |2l|, this one does too, and we have 
shown hidden non-locality. 

We have also tested the new inequality on the 2-qubit 
Werner state l23l: 



Pp =PP[4.) + 



(30) 



where \ip) is the maximally entangled state, as in eqn. 
(|20|l . This state is interesting since despite being entan- 
gled for p > |, Werner gave an explicit LHV model for 
all Von-Neumann measurements for p < i. Ichsh is 
violated for p > , leaving a region | < p < where 
there may or not be model. If there is not a model, it 
must be that some Bell inequality 7iiri222 is violated. We 
find that -Z3322 gives a violation only for p > | , suggesting 
that such a model exists. 

Another important feature of this inequality is that 
non-locality can be shared between qubits. Imagine that 
we have 3 qubits, A, B and C, and we ask whether one 
can simultaneously give non-local correlations with the 
second (summing over the third particle's outcomes), and 
non-local correlations with the third (summing over the 
second particle's outcomes). For CHSH non- locality, the 
answer is "no" ;22|: non-locality is monogamous. We can 
violate the inequality between parties A and B, or A and 
C, or have both pairs give the lhv maximum, but never 
violate both at the same time. 



The non-locality shown by the new inequality can be 
shared. Take the three qubit state 



ABC 



1-f 



(|110) 



ABC 



|101) 



ABC' 



|V) = m|ooo) 

(31) 

with /1 = 0.852. Qubits B and C are symmetric, and 
qubits A and B violate ^3322 giving a value 0.0041. The 
measurements arc defined by the azimuthal and polar 
angles: 



A x = (a, 2tt-/3), 

A 2 = (a,n-/3), 

A 3 = 

B X = ( 7 ,7T + (5), 

B3 = (|,/3), 



(32) 
(33) 
(34) 

(35) 
(36) 

(37) 



where a = 2.8252, (3 = 0.1931, 5 = 0.0804 and 7 = 
2.5445. 

Inspired by this case, we have found a generalization 
of this inequality to the mm-22 case. For 4422 it looks as 
follows: 



^4422 — 
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-1 








o\ 


-3 
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1 


1 


1 


-2 


1 


1 


1 


-1 


-1 


1 


1 


-1 





u 


1 


-1 





0/ 



(38) 



The past-determined correlations are always < 0. The 
generalisation to toto22 should now be clear. The main 
part of the matrix has entries 1 in every position from 
the top left corner to the backwards diagonal. There is 
then one backwards off-diagonal line of —1 entries, and 
afterwards 0's complete the matrix. We then subtract 

Er =1 K-i)p(B i )+p(A 1 ). 

We shall prove the lhv maximum by induction. Start- 
ing from the lower left corner, we can see that Immii con- 
tains all the inequalities from the same family with less 
measurements. If we just take measurements A\ and B m 
(ignoring the other measurements by setting their out- 
comes to 1), we have a positive probability face. Adding 
measurements Ai and B m ^\ gives the CHSH inequality. 
Adding A 3 and B m _ 2 gives us 13322- Let us assume that 
we have proved -f( m _i)(m-i)22 < 0- To get a value larger 
than for l m mTii we must total at least +1 in the terms 
which were not present in the previous inequality. We 
can only do this by setting A4 = OVi, and B\ = 0. This 
gives us +1 in the new terms. But now we have a — 1 
from P(Ai). Whatever we put for the values of Bj for 
j = 2. .rn, each row j contributes exactly to the to- 
tal, giving us I mm 22 — 0, and proving that this is the 
maximum. 

Unfortunately, we do not know if I mm 22 is a face for 
all m. We have found computationally that it is indeed a 
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face for m < 7, and suspect that this will generalise. We 
also do not know if any of the other inequalities mm22 
are relevant after one already has the 3322 inequality. An 
analytic problem here is that we have no simple criterion 
to say which states definitely do not violate the 3322 
inequality. Even computationally we have not yet found 
an example of a state which would violate one of the 
inequalities mm22 without violating 13322- Whilst one 
would expect to find new, inequivalent faces at every m, 
it is not clear that they will all be relevant, particularly 
if we fix ourselves to a certain quantum system size, like 
2-qubits. 

For the case 3422, we find 12480 faces, which include 
three new inequalities, along with 13322, the CHSH and 
positive probability faces. As for our I mm ,22 inequalities, 
we do not know a good way to discover whether these new 
inequalities are relevant, or to uncover other interesting 
features they may possess. The new inequalities are in 
Appendix A. 

We have not gone beyond 3422 in a complete way at 
present since our software takes too long to run on our 
PC. We are able to produce a subset of the faces, but 
have not investigated this direction. 

We can also look at inequalities with more measure- 
ment outcomes. For 2223 there are no new types of face. 
There are only the CHSH, and positive probability faces. 
To use the CHSH inequality (which is defined for two out- 
come measurements) for three outcome measurements, 
we map the three outcomes into two effective outcomes 
by putting two of the original outcomes together. We 
then put the effective outcomes into the CHSH inequal- 
ity. This can be done in 3 different ways for each of 
particle B's measurements. Since there are 8 versions of 
the CHSH inequality for 2222, this gives 8 * 3 ms = 72 
faces. There are mAmsnAnB — 24 positive probability 
faces, making 96 faces in total. 

For 2224 we again find no new types of faces. There are 
32 positive probability faces, and 392 CHSH faces. Note 
that there are two different ways to put together the 4 
outcomes: we can group 3 of them together against the 
4th, or put them in 2 groups of 2. There are 4 ways to do 
the first, and 3 to do the second, giving 8*(4+3) mB = 392 
faces. 

For 2225 and 2226 we have computed a list of all the 
faces, but have not been able to sort them. The number 
of faces, 1840 and 7736, is that which one predicts as- 
suming there are no new faces. Therefore we conjecture 
that there are no faces beyond the CHSH for the case 
222n. 

For the case 2233, there is only one new type of in- 
equality, which is already found in [T2I |24| . This can be 
written as 
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(39) 



The columns of the correlation part of the matrix cor- 
respond to A x = 0, Ax = 1, A 2 = and A 2 = 1. The 
rows are in the same order, for particle B. Thus the first 
entry is P(Ax = 0, Bi = 0). For lhv models I 2 233 < 0. 

The total number of faces for the case 3322 is 1116, of 
which 36 are positive probability, 8 * 3 m ^3 ms — 648 are 
CHSH, and 432 are / 2 233- 

For states of 2 qutrits of the form 



Pp =VP\i,) + (1-Wg, 



(40) 



where 



rk(|o,o> 



li,D 



^vi"'"/ -T- I-, -v t |2,2)), J 2 233 is violated 
by states with more noise (a smaller p) than the CHSH 
inequality. Thus it is relevant. On the other hand, we 
can recover the CHSH inequality from this one by using 
the outcomes 1 and 2 for measurements A\ and A 2 , and 
outcomes and 2 for measurements B\ and B 2 . Once we 
have this new inequality, the CHSH is no longer relevant. 

This inequality has been generalised to 22nn 01 • The 
generalised inequalities are known to be faces for all n, 
and to be the only faces which exist of a certain form (2j| . 
In our present notation, they look simpler than they did 
in the original paper, for example 
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(41) 



The lhv maximum is 0. To see this, note that to get more, 
the local terms —P{B\ ^ n) and —P(Bi ^ n) force us 
to try to get +1 from the three pairs of measurements 
(Ai,Bi), (Ai,B 2 ), and (A 2l B 1 ). But if we do this, we 
are forced to get a —1 from (^2,-82), leaving us with a 
total of 0. 

The generalisation of the inequality to more outcomes 
is as one would guess. 

For 2234 we have computed all the faces, but not sorted 
them. The total number of faces, 19128, matches the 
number one expects assuming there are no new inequal- 
ities. Beyond this our program would take too long to 
compute the complete solution. 

What about inequalities combining more measure- 
ments and more outcomes? Garg and Mermin pr| 
have evidence which suggests that such inequalities ex- 
ist. They found a quantum state and measurements for 
the case 3333 for which the results satisfy all 2222 and 
2233 inequalities, but are nevertheless non-local. Wc 
have found a family of inequalities for the case mmnn, 
which is a generalisation of the inequalities for 22nn and 
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mm22. The first member is 



where X 



J3333 = 
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(42) 



and C 



is a row (or column) in which every entry is C. 



^3333 < for past-generated correlations. X and Y 
are the same matrices which appear in 12233, an d the 
arrangement of X and Y is similar to the arrangement of 
the elements of the matrix in 13322- Z is a new matrix we 
put in by hand, because the more natural matrix full of 
O's did not give us a face. The local probabilities which 
we subtract are a natural generalisation of the terms from 
J2233 and -?3322- 

To generate the complete mmnn family we first gen- 
eralise the number of measurements, then the number of 
outcomes. Adding one more measurement gives 



'4433 
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The generalisation to mm33 follows a similar pattern to 
that for mm22. The matrix consists of entries X for all 
the elements from the top-left corner to just before the 
main backwards diagonal. The main backwards diagonal 
has entries Y, and the next backwards diagonal has en- 
tries —Y. The lower right corner is filled by entries —Z. 
We then subtract P(A 1 £ 2) + ^J 1 {m B -j)P(B j ^ 2). 

To generalise to more outcomes, one only has to change 
the matrices X, Y, and Z. X and Y change exactly as 
they do in the family 22nn. Z grows in a slightly odd 
looking manner: 



Z = 



/oooi\ 

11 
111 
\0 0/ 



(44) 



which is the same as Y but for the last row. 

Immnn < for lhv theories. To prove this, we combine 
the proofs of I mm i2 < and I 22nn < 0. Starting at the 
bottom left corner, and successively adding pairs of mea- 
surements, we see that I mmnn contains all the inequali- 
ties Im'm'nn, with m! < m. For m = l, the inequality 
is trivial. For m = 2, the inequality is I 22nn , which we 
have already proved. For m = 3, to get more than we 
need to get a positive contribution from the terms added 
on after m = 2. Thus we must get a +1 from all the 
combinations (A\,B\), (A 2 ,Bi), (A 3 ,Bi). But now we 



have —1 from A\. To get a total of more than 0, we need 
to find a row Bk where the contribution from that row is 
+1. To see that this is impossible, first look at the row 
B 2 . We want to pick up +1 from (Ai,B 2 ) and (A 2 ,B 2 ) 
without picking up a —1 from (A 3 ,B 2 ). But the +l's in 
(^2,-62), (A 2 , B\) and (A 3 ,Bi) force the —1 to occur, 
making the maximum of the row 0. Otherwise look at 
row B 3 . Here setting B 3 = n — 2 is no use, since we get 
a +1 and a— 1. B3 = n — 1 is clearly useless, also giving 
us 0. For B 3 < n — 2, the matrix — Z looks like —Y, 
and so will always give —1 for (A 3 , B 3 ) when we get +1 
for (Ai,B 3 ), (A\,Bi) and {A 3 ,Bi). So we have proven 
^33nn < 0. For more measurements, a similar argument 
leads to a proof by induction. 

We know computationally that these inequalities are 
faces for m = 2, n < 7, for m = 3, n < 6, for m — 
4, n < 4, and for m = 5, n = 3. We suspect this is 
true for all m and n. As was the case for our family of 
mm22 inequalities, we do not know if these new mmnn 
inequalities are relevant. 

In summary we have found that for small numbers of 
measurements and outcomes there are very few inequiv- 
alent Bell inequalities. For the case 2m22 there is only 
the CHSH inequality. We believe the same to be true for 
the case 222n. For the cases 2233 and 3322 there is the 
CHSH inequality, but only one other inequality in each 
case. The new inequalities are relevant: they are violated 
by states which do not violate the CHSH inequality. The 
3322 case shows that the CHSH inequality is not the only 
useful one for 2-qubits: more measurements really help! 
We have also discovered a family of inequalities for the 
case mmnn, but do not know if any of these inequalities 
are relevant. 

Having found that there are remarkably few relevant 
incquivalent faces, we have more time to study closely 
the ones which do exist. Will some of them help us to 
perform an experiment definitively ruling out lhv correla- 
tions in the lab? Is there a close connection between the 
new inequalities and a particular quantum cryptography 
protocol? How does this compare with different types of 
quantum correlations? All these possibilities would be 
interesting, but here we have a more fundamental mes- 
sage. The set of past-generated correlations is not as 
complicated as we thought. 

Thanks to: A. Acin, N. Brunncr, R. Gill, N. D. Mer- 
min, I. Pitowsky, V. Scarani and A. Shimony for helpful 
comments and suggestions, and S. Fasel for software sup- 
port. 

We acknowledge funding by the Swiss NCCR, "Quan- 
tum Photonics" and the European 1ST project RESQ. 

Note Added: C. Sliwa has independently found some 
of the results contained in this paper. Firstly that for the 
case 2n22 there are no inequalities beyond the CHSH. 
Secondly that for the case 3322 there is only a single new 
inequality (this is, like our result, an exact computational 
result). He has also investigated the three party case 
222222. 
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APPENDIX A: NEW 3422 INEQUALITIES 

There are 12480 faces for the case 3422, which include 
3 new faces. There are 3*4*2*2 = 48 positive probability 

faces, 8 (^j (^j = 144 CHSH faces ' and 576 ( 3 ) = 

2304 ^3322 faces. Then the three new inequalities. There 
are 2304 versions of 
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which has a lhv maximum of 2. There are 3027 versions 
of 
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with a lhv maximum of 1. Finally 4608 versions of 
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with a lhv maximum of 2. 
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